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Abstract. Adaptive, locally refined and locally adjusted meshes are pre- 
ferred over uniform meshes for capturing singular or localised solutions. 
Roughly speaking, for a given degree of freedom a solution associated 
with adaptive, locally refined and locally adjusted meshes is more accu- 
rate than the solution given by uniform meshes. In this work, we answer 
the question which meshes are better conditioned. We found, for approx- 
imately same degree of freedom (same size of matrix), it is easier to solve 
a system of equations associated with an adaptive mesh. 



1 Introduction 

Uniform, locally adjusted, adaptive and locally refined meshes are shown in 
Figures |3 and El respectively. Here, for each mesh the number of cells (or 
degree of freedom) are approximately 1024 (2 5 x 2 5 ). Let us consider the steady 
state pressure equation of a single phase flow in a porous medium fl 

— div (K gradp) = / in fl and p(x,y)=p D on dfl^ ■ (1) 

Here, fl is a polyhedral domain in R 2 , the source function / is assumed to 
be in L 2 (fl) and the diagonal tensor coefficient K(x,y) is positive definite and 
piecewise constant. K (permeability) is allowed to be discontinuous in space. We 
are discretizating the equation on the meshes (see Figures |3 01 \5\ and EJ by 
the method of Finite Volumes Ill.'il7l8| . For discretization of the problem (JIJ on 
uniform and localised meshes (see the Figures|21and01) , we refer to the References 
|1I7I8| . Discretization of the equation (JIJ on adaptive and locally refined meshes 
is given in the following References |3I5| . Finite Volume discretization of the 
problem Qona mesh results in a matrix system Aph = b. Here, A is symmetric 
positive definite matrix associated with a mesh. 

Let us define a problem to be solved on the four meshes. Let the domain be 
fl = [—1, 1] x [—1, 1] (see Figure^. It is divided into four sub-domains according 
to the permeability K (see the Figures ^ and 2J. The permeability K is a positive 
constant in each of the sub-domains and is discontinuous across the surfaces of 
sub-domains. Let the permeability in the sub-domain fli be K^. Assuming that 
Ki = K3 = R and K2 = K4 = 1.0. Ki, K2, K3 and K4 refers to the permeabilities 



in the subdomains fix, J?2 ; ^3 and Q4, respectively. The parameter R is given 
below. Let the exact solution in the polar form be [3] 

p(r,9)=r^r,(9) , (2) 

where the parameter 7 denotes the singularity in the solution [5] and it depends 
on the permeability distribution in the domain (see Figure^for the permeability 
for the singularity 7 = 0.1). r](d) is given as 

'cos[(tt/2 - cr) 7 ] cos[(0 - vr/2 + j o) 7 ] , 9 e [0, tt/2] , 
/m = 1 cos(/?7) cos[(0 - 7T + o-) 7 ] , 9 e [tt/2, tt] , 

C0S((77) COS[(0 — 7T — /?)7] , 6* G [7T, 37T/2] , 

k cos[(tt/2 - p)7] cos[(6> - 3vr/2 - cr) 7 ] , (9 e [3tt/2, 2tt] . 

It can be shown that solution p (given by equation (0)) barely belongs in the 
fractional Sobolcv space H 1+K (i?) with n < 7 (cf. [5]). 
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Fig. 1. Domain. Fig. 2. Permeability distribution. 



For the singularity 7 = 0.1, the parameters are 

R w 161.4476 , p w 0.7854 and a w -14.9225 . 

We solve the problem on the four meshes. The exact solution is given by 
the equation @. We enforce the solution inside the domain by the Dirichlet 
boundary condition and the source term. For solving discrete system of equations 
formed on the meshes, we use the Conjugate Gradient (CG) solver (see 0]). Table 
n presents eigenvalues and condition numbers of the matrix systems associated 
with the different meshes. Note that in this table, the largest eigenvalue on all 
four meshes is approximately same. However, the smallest eigenvalue associated 
with the adaptive mesh is greater than the smallest eigenvalues associated with 
other three meshes. When solving the Symmetric Positive Definite (SPD) linear 
system Ap/, = b with the CG, the smallest eigenvalues of the matrix slowes 
down the convergence (cf. T). Several techniques have been proposed in the 
literature to remove bad effect of the smallest eigenvalue (see [2141 and references 
therein]). Convergence of the CG solver for these the four systems are shown in 
the Figure [7| It is clear from the Table ^ and the Figure [7| that it is easier to 
solve a matrix system associated with an adaptive mesh than to solve systems 
associated with uniform, localised and locally refined meshes. 



Fig. 3. Uniform mesh. 



Fig. 4. Localised mesh. 




Fig. 5. Adaptive mesh. Fig. 6. Locally refined mesh. 




Iterations 

Fig. 7. Convergence curves for the matrix system formed on the different meshes. 



Table 1. Eigenvalues and condition numbers of different matrix systems. 



Mesh 


Smallest eig. 


Largest eig. 


Cond. Num. 


Adaptive 


4.15 x 10" 1 


1.28 x 10 3 


3.10 x 10 3 


Localised 


5.50 x 1CT 2 


0.78 x 10 3 


1.42 x 10 4 


Uniform 


7.62 x f(T 2 


1.28 x 10 3 


1.69 x 10 4 


Locally Refined 


3.94 x f (T 2 


1.28 x 10 3 


3.25 x 10 4 



2 Conclusions 

We have shown that it is easier to solve a matrix system associated with an 
adaptive mesh than solving systems associated with uniform, localised and lo- 
cally refined meshes. The adaptive mesh is generated by equal distribution of 
the fluxes over all the cells in the mesh. Why do equal distribution of fluxes is 
create meshes which are better conditioned ? Or, why do equal distribution of 
fluxes remove bad effect of small eigenvalue ? Answers to these questions can 
help in designing new preconditioners or improving existing preconditioners. 
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Preface 



This textbook is intended for use by students of physics, physical chemistry, 
and theoretical chemistry. The reader is presumed to have a basic knowledge 
of atomic and quantum physics at the level provided, for example, by the first 
few chapters in our book The Physics of Atoms and Quanta. The student of 
physics will find here material which should be included in the basic education 
of every physicist. This book should furthermore allow students to acquire an 
appreciation of the breadth and variety within the field of molecular physics and 
its future as a fascinating area of research. 

For the student of chemistry, the concepts introduced in this book will provide 
a theoretical framework for that entire field of study. With the help of these con- 
cepts, it is at least in principle possible to reduce the enormous body of empirical 
chemical knowledge to a few basic principles: those of quantum mechanics. In 
addition, modern physical methods whose fundamentals are introduced here are 
becoming increasingly important in chemistry and now represent indispensable 
tools for the chemist. As examples, we might mention the structural analysis of 
complex organic compounds, spectroscopic investigation of very rapid reaction 
processes or, as a practical application, the remote detection of pollutants in the 
air. 
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1 Fixed-Period Problems: The Sublinear Case 

With this chapter, the preliminaries are over, and we begin the search for periodic 
solutions to Hamiltonian systems. All this will be done in the convex case; that 
is, we shall study the boundary-value problem 

x = JH'(t,x) 
x(0) = x(T) 

with H(t, •) a convex function of x, going to +00 when ||x|| — ► 00. 
1.1 Autonomous Systems 

In this section, we will consider the case when the Hamiltonian H(x) is au- 
tonomous. For the sake of simplicity, we shall also assume that it is C x . 

We shall first consider the question of nontriviality, within the general frame- 
work of (^4oo, i?oo)-subquadratic Hamiltonians. In the second subsection, we shall 
look into the special case when H is (0, 6 oc )-subquadratic, and we shall try to 
derive additional information. 

The General Case: Nontriviality. We assume that H is {A^, _B oc )-sub- 
quadratic at infinity, for some constant symmetric matrices A^ and B^, with 
Boo — Aoo positive definite. Set: 



7 : = smallest eigenvalue of — A, 



00 



(1) 



A : = largest negative eigenvalue of J— + A, 



dt 



OO 



(2) 



Theorem ?? tells us that if A + 7 < 0, the boundary- value problem: 

x = JH'{x) 
x(0) = x{T) 



(3) 



has at least one solution x, which is found by minimizing the dual action func- 
tional: 

rT 



±(A?u,u)+N*{-u) 



dt 



(4) 



on the range of A, which is a subspace R(A)\ with finite codimension. Here 



N{x) := H{x) - - {A^x^x) 



is a convex function, and 

N(x) < i ((Boo -Aoo)x,x) +c \fx . 
Proposition 1. Assume H'(0) = and H(0) = 0. Set: 



S := liminf 2N(x) \\x 

7/7 < —A < S, the solution u is non-zero: 

x{t) ± Vt . 



1-2 



(5) 



(6) 



(7) 



(8) 



Proof. Condition (??) means that, for every 8' > 6, there is some e > such 
that 



6'.. ,.2 



\\x\\ < e => N(x) < — \\x 



(9) 



It is an exercise in convex analysis, into which we shall not go, to show that 
this implies that there is an r\ > such that 



f\\x\\< V ^N*(y)<—\\y 



(10) 



Fig. 1. This is the caption of the figure displaying a white eagle and a white horse on 
a snow field 



3 



Since u\ is a smooth function, we will have H/imiH^ < i] for h small enough, 
and inequality (??) will hold, yielding thereby: 



ip(hut) < \\ \\ Ul \\l + ^-^ ||ui|| 2 



(11) 



If we choose 5' close enough to S, the quantity ( j + jr) will be negative, and 
we end up with 

ij){hu{) < for h ^ small . (12) 

On the other hand, we check directly that tp(0) = 0. This shows that cannot 
be a minimizer of tjj, not even a local one. So u ^ and u ^ A~ 1 (0) = 0. □ 

Corollary 1. Assume H is C 2 and (doo, &oo) -subquadratic at infinity. Let £i, 
...,£n be the equilibria, that is, the solutions of H'(£) = 0. Denote by uik the 
smallest eigenvalue of H" (£fc), and set: 



uo := Min {u)\, . . . , w*;} 



T 



(13) 
(14) 



then minimization of tp yields a non-constant T -periodic solution x. 



We recall once more that by the integer part E[a] of a G H, we mean the 
a e 7Z such that a < a < a + 1. For instance, if we take = 0, Corollary 2 
tells us that x exists and is non-constant provided that: 



2^ 6 ~ <1< 2^ 



or 



T G 



2tt 2tt 



(15) 
(16) 



Proof. The spectrum of A is 1^.2? + a^. The largest negative eigenvalue A is 
given by + a^, where 



— A: + aoo < < — (k„ + 1) + a c 



Hence: 

The condition 7 < —A < S now becomes: 

2tt 



aoo < 



k a - < lu - a a 



which is precisely condition (??). 



(17) 
(18) 

(19) 
□ 
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Lemma 1. Assume that H is C 2 on R 2 ™\{0} and that H"(x) is non-degenerate 
for any x ^ 0. Then any local minimizer x of ip has minimal period T. 

Proof. We know that x, or x + £ for some constant £ G H 2 ™, is a T-periodic 
solution of the Hamiltonian system: 

x = JH'(x) . (20) 

There is no loss of generality in taking £ = 0. So ip(x) > i[)(x) for all x in 
some neighbourhood of x in IT 1 ' 2 (TR/T2Z; IR 2 ™). 

But this index is precisely the index ix{x) of the T-periodic solution x over 



the interval (0,T), as defined in Sect. 2.6. So 

ir(x) - . (21) 
Now if x has a lower period, T/k say, we would have, by Corollary 31: 

i T (x) = i kT/k (x) > ki T/k (x) + k-l>k-l>l . (22) 
This would contradict (??), and thus cannot happen. □ 



Notes and Comments. The results in this section are a refined version of [?] ; the 
minimality result of Proposition 14 was the first of its kind. 

To understand the nontriviality conditions, such as the one in formula (??), 
one may think of a one-parameter family xt, T G (27rw _1 , 2irb^) of periodic 
solutions, xt(0) = %t(T), with xt going away to infinity when T — > 2irui~ 1 , 
which is the period of the linearized system at 0. 

Table 1. This is the example table taken out of The T^Xbook, p. 246 



Year 


World population 


8000 B.C. 


5,000,000 


50 A.D. 


200,000,000 


1650 A.D. 


500,000,000 


1945 A.D. 


2,300,000,000 


1980 A.D. 


4,400,000,000 



Theorem 1 (Ghoussoub-Preiss). Assume H(t,x) is {Q 1 e)-subquadratic at 



infinity for all e > 0, and T-periodic in t 

H(t, •) is convex Vi (23) 

H(-,x) is T-pcriodic Vx (24) 

H(t,x) > n(\\x\\) with n(s)s^ 1 — > oo as s — > oo (25) 
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Ve>0, 3c : H(t,x) < ^\\xf + c . (26) 

Assume also that H is C 2 , and H"(t, x) is positive definite everywhere. Then 
there is a sequence Xk, k <E IN, of kT -periodic solutions of the system 

x = JH'(t,x) (27) 

such that, for every k G IN, there is some p a € IN with: 

P>Po=>x pk ^x k . (28) 

□ 

Example 1 (External forcing ). Consider the system: 

x = JH'{x) + f(t) (29) 

where the Hamiltonian H is (0, &oo)-subquadratic, and the forcing term is a 
distribution on the circle: 

f=j t F + f with FeL 2 (M/T2Z; IR 2 ™) , (30) 
where f a := T _1 f(t)dt. For instance, 

/(*) = E s *z > ( 31 ) 

feeiN 

where S k is the Dirac mass at t = k and £ G IR 2 " is a constant, fits the pre- 
scription. This means that the system x = JH'(x) is being excited by a series 
of identical shocks at interval T. 

Definition 1. Let A^it) and i?oo(£) be symmetric operators in JR 2n , depending 
continuously on t G [0, T], such that ^(t) < B^t) for all t. 

A Borelian function H : [0, T] x JR 2n — > IR is called (A^, B^-subquadratic 



at infinity if there exists a function N(t, x) such that: 

H(t,x) = ^{A 00 {t)x,x)+N(t,x) (32) 

Vi , N(t, x) is convex with respect to x (33) 

N(t,x) > n(\\x\\) with n(s)s^ 1 -> +oo as s -> +oo (34) 

BcelR: H(t,x)<^(B 00 (t)x,x)+c Vx . (35) 



If A^t) = aool and B^t) = b^I , with < G 1R ; we shall say that 
H is (aocb^-subquadratic at infinity. As an example, the function \\x\\ a , with 
1 < a < 2, is (0, e)-subquadratic at infinity for every e > 0. Similarly, the 
Hamiltonian 

H{t,x) = ^k\\k\\ 2 + \\x\\ a (36) 
is (k, k + e)-subquadratic for every e > 0. Note that, if k < 0, it is not convex. 
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Notes and Comments. The first results on subharmonics were obtained by Ra- 
binowitz in [?] , who showed the existence of infinitely many subharmonics both 
in the subquadratic and superquadratic case, with suitable growth conditions 
on H'. Again the duality approach enabled Clarke and Ekeland in [?] to treat 
the same problem in the convex-subquadratic case, with growth conditions on 
H only. 

Recently, Michalek and Tarantello (see [?] and [?]) have obtained lower bound 
on the number of subharmonics of period kT, based on symmetry considerations 
and on pinching estimates, as in Sect. 5.2 of this article. 
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1 Fixed-Period Problems: The Sublinear Case 

With this chapter, the preliminaries are over, and we begin the search for periodic 
solutions to Hamiltonian systems. All this will be done in the convex case; that 
is, we shall study the boundary-value problem 

x = JH'(t,x) 
x(0) = x(T) 

with H(t, •) a convex function of x, going to +oo when ||x|| — ► oo. 
1.1 Autonomous Systems 

In this section, we will consider the case when the Hamiltonian H(x) is au- 
tonomous. For the sake of simplicity, we shall also assume that it is C . 

We shall first consider the question of nontriviality, within the general frame- 
work of (Am, -Boo)-subquadratic Hamiltonians. In the second subsection, we shall 
look into the special case when H is (0, fe^-subquadratic, and we shall try to 
derive additional information. 

The General Case: Nontriviality. We assume that H is {A^, B^-sub- 
quadratic at infinity, for some constant symmetric matrices A^ and B^, with 
Bqc — Aoo positive definite. Set: 



7 : = smallest eigenvalue of B x — A, 



oc 



(1) 



A : = largest negative eigenvalue of J— + A, 



dt 



(2) 



Theorem 21 tells us that if A + 7 < 0, the boundary-value problem: 



x = JH'(x) 
x(0) = x(T) 



(3) 



has at least one solution x, which is found by minimizing the dual action func- 
tional: 

" ~1 



ip(u) 



-(A^u,u)+N*(-u) 



dt 



(4) 



on the range of A, which is a subspace R(A) 2 L with finite codimension. Here 



N(x) :=H(x)--(A oc x,x) 



is a convex function, and 



(5) 



N(x) <-((B 00 -A 00 )x,x)+c Vx 

Proposition 1. Assume H'(0) = and H(0) = 0. Set: 

5 := lim inf 2N(x) \\x\\~ 2 . 

If 7 < —A < d, the solution u is non-zero: 

x{t) ^0 Vt . 



(6) 



(7) 



(8) 



Proof. Condition (??) means that, for every 8' > S, there is some e > such 
that 

\\x\\<e^N(x)< S -\\x\\ 2 . (9) 

It is an exercise in convex analysis, into which we shall not go, to show that 
this implies that there is an r\ > such that 



fM<V^N*(y)<±\\yf 



(10) 



Fig. 1. This is the caption of the figure displaying a white eagle and a white horse on 
a snow field 

Since u\ is a smooth function, we will have H/iuiH^ < f] for h small enough, 
and inequality (??) will hold, yielding thereby: 

^i)<yi|hii||^ + y^|| U i|| 2 • (11) 
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If we choose 6' close enough to 6, the quantity ( j + jr) will be negative, and 
we end up with 

ipihm) < for h ^ small . (12) 

On the other hand, we check directly that ip(0) = 0. This shows that cannot 
be a minimizer of ■0, not even a local one. Som^O and u ^ A~ 1 (0) = 0. □ 

Corollary 1. Assume H is C 2 and {a^, boo) -subquadratic at infinity. Let £i, 
. . . , £jv be the equilibria, that is, the solutions of H' '(£) = 0. Denote by uik the 
smallest eigenvalue of H" (£k), an d se t : 



:= Min {ui, . . . ,Lu k } . 



If: 



2tt 



2tt 



T 



(13) 
(14) 



then minimization of tp yields a non-constant T -periodic solution x. 



We recall once more that by the integer part E[a] of a G IR, we mean the 
a e 2Z, such that a < a < a + 1. For instance, if we take Ooc = 0, Corollary 2 
tells us that x exists and is non-constant provided that: 



or 



T , T 
— boo < 1 < ^~ 
2tt 2tt 



/ 2tt 2vr 

Te — r- 



(15) 



(16) 



Proof. The spectrum of yl is + a^. The largest negative eigenvalue A is 

given by ^fc Q + a^, where 



2tT 27T., 

— Ko + doo < < — \ k o + 1) + a c 



Hence: 



k a = E 



T 



The condition 7 < —A < 5 now becomes: 



2tt, 



which is precisely condition (??). 



(17) 
(18) 

(19) 
□ 



Lemma 1. Assume that H is C 2 on TR 2n \{0} and that H"(x) is non-degenerate 
for any x ^ 0. Then any local minimizer x of ip has minimal period T. 
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Proof. We know that x, or x + £ for some constant £ e H n , is a T-periodic 
solution of the Hamiltonian system: 

x = JH'(x) . (20) 

There is no loss of generality in taking £ = 0. So ip(x) > ip(x) for all x in 
some neighbourhood of x in W 1 ' 2 (M/TZZ; H 2 ™). 

But this index is precisely the index ir{x) of the T-periodic solution x over 
the interval (0, T), as defined in Sect. 2.6. So 

ir(x) = . (21) 

Now if x has a lower period, T/k say, we would have, by Corollary 31: 

*t(x) = i kT /k{x) > ki T/k (x) + k - 1 > k - 1 > 1 . (22) 

This would contradict (??), and thus cannot happen. □ 

Notes and Comments. The results in this section are a refined version of ?; the 
minimality result of Proposition 14 was the first of its kind. 

To understand the nontriviality conditions, such as the one in formula (??), 
one may think of a one-parameter family xt, T e (27ra; -1 , 27r&^ 1 ) of periodic 
solutions, Xt(0) = %t(T), with xt going away to infinity when T — > 2iruj~ 1 , 
which is the period of the linearized system at 0. 



Table 1. This is the example table taken out of The TgXbook, p. 246 



Year 


World population 


8000 B.C. 


5,000,000 


50 A.D. 


200,000,000 


1650 A.D. 


500,000,000 


1945 A.D. 


2,300,000,000 


1980 A.D. 


4,400,000,000 



Theorem 1 (Ghoussoub-Preiss). Assume H(t,x) is (0, s)-subquadratic at 
infinity for all e > 0, and T-periodic in t 

H(t, •) is convex Vt (23) 

H(-,x) is T-pcriodic Vx (24) 
H(t,x) > n(\\x\\) with n(s)s^ 1 — > co as s — > co (25) 

Ve>0, 3c : H(t,x) <^\\xf + c . (26) 
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Assume also that H is C 2 , and H"(t, x) is positive definite everywhere. Then 
there is a sequence Xk, fc £ IN, of kT -periodic solutions of the system 

x = JH'(t,x) (27) 

such that, for every k £ IN, there is some p £ IN with: 

P> Po => x pk ^ x k . (28) 

□ 

Example 1 (External forcing ). Consider the system: 

x = JH'(x) + f(t) (29) 

where the Hamiltonian H is (0, &oo)-subquadratic, and the forcing term is a 
distribution on the circle: 

f=jF + f Q with F £ L 2 (M/T2Z; JR 2n ) , (30) 
where f := T _1 f(t)dt. For instance, 

/(*) = E s *z - ( 31 ) 

feeiN 

where 5k is the Dirac mass at t = k and £ £ M 2 ™ is a constant, fits the pre- 
scription. This means that the system x = JH'(x) is being excited by a series 
of identical shocks at interval T. 

Definition 1. Let A^t) and B^t) be symmetric operators in IR 2 ™, depending 
continuously on t £ [0, T], such that ^(t) < B^t) for all t. 

A Borelian function H : [0, T] x IR 2 " — > IR is called (A^, B^-subquadratic 



at infinity if there exists a function N(t, x) such that: 

H(t,x) = ^{A 00 (t)x,x)+N(t,x) (32) 

Vt , N(t, x) is convex with respect to x (33) 

N(t,x) > n(\\x\\) with n(s)s^ 1 -> +oo as s -> +oo (34) 

3c£H: H(t,x) < ^ (Boo(t)a;,a;) +c Vx . (35) 



If Aoo(t) = a^I and B^t) = b^I , with < b^ £ IR, we shall say that 
H is (aaojbooj-subquadratic at infinity. As an example, the function \\x\\ , with 
1 < a < 2, is (0, e)-subquadratic at infinity for every e > 0. Similarly, the 
Hamiltonian 

H(t,x) = h\\kf + \\x\\ a (36) 
is (k, k + e)-subquadratic for every e > 0. Note that, if k < 0, it is not convex. 
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Notes and Comments. The first results on subharmonics were obtained by Ra- 
binowitz in ?, who showed the existence of infinitely many subharmonics both 
in the subquadratic and superquadratic case, with suitable growth conditions 
on H'. Again the duality approach enabled Clarke and Ekeland in ? to treat 
the same problem in the convex-subquadratic case, with growth conditions on 
H only. 

Recently, Michalek and Tarantcllo (see Michalek, R., Tarantello, G. ? and 
Tarantcllo, G. ?) have obtained lower bound on the number of subharmonics of 
period kT, based on symmetry considerations and on pinching estimates, as in 
Sect. 5.2 of this article. 
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